ADA053521 


NORTHWESTERN- UNIVERSITY  | 

.The  TechnologicJil  Institute  , | 

Department  of  Engineering  i 

Evanston,  Illinois  60201'  • | 


THE  STRESS  FIELD  CREATED  BY  A 
CIRCL'L.\R  SLIDING  CONTACT 
ON  TRANS\^ER5ELY  ISOTROPIC  SFHERES 


A THESIS 

SUBMITTED  TO  THE  GRADUATE  SCHOOL 
IN  PARTIAL  FULJTLLMEKT  OF  THE  REQUIREMENTS 


for  the  degree 


MASTER  OF  SCIENCE 
Field  of  Civil  Engineeririg 


Dale  Bruce  Mowrj- 


E^'nn^■ton,  Illinois 


DISTRIBUTION  STJffEMENT  M 

Approved  for  public  xeleoeej 
DiiizibutkMi  (AoJiodted 


June  197S 


D D 

SCnnnET 

WAY  4 19TS 


NOmWESTERN  UNIVERSIT7 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

1 


SUBMITTED  TO  THE  GRADUATE  SCHOOL 
IN  PARTIAL  FULFILIMENT  OP  THE  REQUIREMENTS 


for  the  degree 

MASTER  OF  SCIENCE 
Field  of  Civil  Engineering 

by 

le^  Bruc^ Mowry^ 
Evanston,  Illinois 


@[ 


Acknowledgments 


I wish  to  express  hqt  sincere  gratitude  to  Dr.  Leon  M.  Keer, 
Professor  of  Civil  Engineering,  Northwestern  University,  tdio  has  given 
generously  of  his  time  and  guidance  in  the  completion  of  this  project, 
and  to  oy  wife,  Nancy,  for  her  loving  encouragement  and  understanding. 


Abstract 


\ 

The  field  equations  for  a transversely  isotropic  half>space  are 
defined  in  tezms  of  potential  functions.  A form  of  the  potential 
functions  satisfying  equilibrium  is  assumed  and  the  boundary  value  prob- 
lem of  a tangentially  loaded  spherical  indentor  on  the  half -space  is  | 

} 

solved.  ' 

Expressions  are  obtained  defining  the  radius  of  no-slip  for  the 
static  case,  and  the  relationship  between  the  horizontal  surface  dis- 
placement under  the  Indentor  and  the  horizontally  applied  force. 

Stresses  for  tangential  loading  are  superposed  with  diose  previously 
obtained  for  normal  loading  of  the  indentor,  and  the  stress  field  is 
defined  in  the  half- space  and  on  the  surface  for  both  static  and  slid- 
ing cases. 

Von  Hlses'  criteria  for  the  sliding  case  is  calculated  and  plotted 
in  the  half-space  and  on  the  surface  for  two  transversely  isotropic 
metals  using  two  separate  coefficients  of  limiting  friction,  and  on  the 
surface  for  the  static  case  for  magnesium  with  a coefficient  of  limiting 
friction  of  .5.^ 
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1.  Introduction 


The  stress  field  created  by  a circular  sliding  contact  on  an  isotop- 
ic half-space  has  been  solved  by  Hamilton  and  Goodman  [1].  Motivated  by 
consideration  of  mechanical  failure  they  examined  constant  lines  of 
von  Mises'  yield  criteria  in  the  half -space  and  on  the  surface.  Dahan 
and  Zarka  [2]  have  recently  solved  the  stress  field  in  a transversely 
isotropic  half-space  in  contact  with  a spherical  Indentor  under  normal 
loading.  They  also  plotted  von  Mises*  criteria  for  several  transversely 
isotropic  metals  to  show  the  effect  of  the  anisotropy  for  the  indentation 
of  an  elastic  half-space.  Both  solutions  give  results  in  the  half-space 
in  a closed  form. 

The  purpose  of  the  present  analysis  is  to  derive  the  stress  ex- 
pressions for  identical  transversely  isotropic  spheres  in  contact 
under  both  normal  and  tangential  loading,  and  to  obtain  expressions  for 
the  stresses  as  the  spheres  slide  relative  to  each  other  as  a result  of 
the  tangential  load.  Von  Mises*  yield  criteria  is  then  examined  to  com- 
pare the  anisotropic  case  to  the  isotropic  case.  iA  [IJ. 
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2.  FoxBulatlon 


(a)  Basic  Equations 

For  a transvarsaly  Isotropic  body  the  stress-strain  relations  are 
given  by  Green  and  Zema  [3]  as 


" ®11  ®12  ®13  ° ° 
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11  13 


c 0 0 

33 


®44  ° 


e I (1) 

yz 


%(c  -c  ) e 
11  12  xy 


The  relations  between  the  strains  and  displacements  are 

3u  3u  au^ 

e - — - , e - — i , e - 

« Bx  yy  3y  « a* 


yz  ' a*  ay/  **  '3*  ax/  xy  \ ay  ax/ 


and  the  stress-strain  relations  become 
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Substitution  of  the  above  stress  equations  Into  the  equilibrium 
equations,  aa  /ax  « 0,  gives  the  following  expressions  for  equilibrium 

Ij  j 

in  terms  of  displacement: 


3x^ 


+ if  ] ■ ® 


+ a/ay  [«'=ii+<=i2>  ^ + <'u«44>  ^ ] * ® 


ax^  ay^^^ 
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The  displacements  are  defined  next  by  potential  functions  as  in  [3]: 


u ■ + 1^  , u-  ■ - 5JL  , u * k 

* 3x  ay  ^ ay  ax  * az 


idiere  k Is  a constant  to  be  determined  later. 

Substituting  these  displacement  expressions  into  the  equilibrium 
equations  (4)  leads  to 
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Clearly  these  equations  are  satisfied  If 


* t'44  * ''<'i3^44>]  ^ - » 
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2.  2 


“u  ‘44  * •“ 
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Equations  (7a)  and  (7b)  can  be  written 


(ijL  + A)  + ‘44  * ‘‘('=U  * W ^ j 
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As  pointed  out  In  [3],  a suitable,  non-zero  0 can  be  found  If  (8)  and 
(9)  are  Identical,  vdilch  leads  to 
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Equation  (10)  generates  a quadratic  equation  for  the  solution  of  v , 


®11®44''  ^®13^^®44  ”*■  ®13^  " ®11®33^  * *^33^44  " ^ 


from  which  It  Is  found  that 
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with  the  functions  0^  and  0^  applicable  to  Vj^  and  v^*  Hence  the  equations 


that  determine  0 and  0.  are 
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Fron  (15)  and  (16)  a form  of  the  potential  functions  Is  assuned  as 
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, ? -§2//v 
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n»0  o 
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To  determine  Y,  equation  (7c)  is  written  in  the  form 
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Then, 
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and  0 has  the  form 
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^ • I I *3„<« 


-|2//v. 
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The  form  of  the  displacement  is 
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with  (17),  (18),  «nd  (21)  representing  the  functions  0^,  0^  and 


0 the  application  of  the  boundary  conditions  determines  A (C) . d . 
3 in  ^ ^i' 


and  finally,  from  the  stress-strain  relations  (1),  the  stress  field. 

Ihe  stress  field  for  the  case  of  two  Identical  transversely 
isotropic  spheres  in  contact  is  obtained  by  first  solving  the 
problem  of  normally  loaded  spheres.  Then  the  case  of  tangential  load 
in  the  positive  x-dlrection  insufficient  for  sliding  to  occur  is  con- 
sidered, where  sliding  refers  to  the  condition  of  total  slip  between 
the  two  spheres.  The  stress  field  obtained  is  then  specialized  to 
give  results  for  the  case  where  sliding  of  the  sphere  occurs.  The 
problem  of  normally  loaded  transversely  isotropic  spheres  as  previously 
noted,  has  been  solved  in  [2]  so  that  the  following  solution  is  for 
tangential  loading.  The  two  solutions  will  then  be  superposed. 

with  the  orientation  of  the  axes  as  shotm  in  figure  2.1,  the  boundary 
conditions  for  the  application  of  a tangential  force  in  the  positive 


x-directton,  P^,  are  (when  appropriate  symmetry  conditions  are  considered 


for  the  two  spheres) 


3 


f 


z 


4 o 

'z 


Figure  2.1  Physical  Orientation  of  the  Problem 
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on  z ■ 0: 


u » 0 

y 


0 < r < c 


u ■ a 

X 


0 < r < c 


a » g(r) 
zx 


c < r < a 


O’  *0 
zx 


a < r < oa 


o *0 
zy 


c < r < SB 


where  0 < r < £ is  the  region  of  no  slip,  c < r < a is  the  region  diat 
experiences  slip,  o is  constant,  and  g(r)  is  a function  of  r to  be 
determined  later. 

Since  the  forms  of  the  potentials  are  In  cylindrical  polar  coordinates. 


the  boundary  conditions  are  put  In  that  coordinate  system  as  follows: 


u > u cos  6 + u sin  @ « o cos  6 
r X y 


u ■ u cos  9 - u sin  9 ■ -a  sin  9 
9 7 * 


0 < r < c 


0 < r < £ 


If  u^  and  u^  are  expressed  In  Fourier  cosine  and  sine  series, 


u > ; u cos  n6  ■ or  cos  9 

r ^ m 

n«»0 


0 < r < £ 


u-  ■ fi  sin  n9  ■ *«  sin  9 
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0 < r < £ 
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From  ayuiMtry  considerations  and  equation  (25),  n » 1 and 


u » u cos  Q 
r r 


u « u sin  9 

s e 


Then  In  the  region  0 < r < c , 


“x  " "r  e - Ug  sin  9 - - fi  ) + %(a^  + d ) cos  29 


u « u sin  6 + u cos  9 - %(u  + 6)  sin  28 
y ^ H r 0 


From  (27a)  and  (23),  with  A - 2a, 


A A . 

U - U - A , 

r 9 


a + fi  - 0 , 
r 8 


0 < r < c 


0 < r < 6 


In  the  same  manner,  with  a • 6 cos  9,  o » a sin  9, 

zr  zr  z6  z9 

and  f(r)  - 2g(r), 


(27a) 
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(28b) 


“ ^^"^zr  ” ^z9^  ^^^zr  ^z9^  ^ " 8(r),  £ < r < a 


yields 
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*1^  z9 
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d + a - 0 , 
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(31b) 
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Since 


on  z «*  0. 
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(32a) 


and  a ■ %(a  + a ) 

zy  zr  z8 

sin  26, 

(32b) 

then  <?  + ^ - “ 0, 

zr  z0 

5 < r < o» 

(33a) 

(T  - CT  - 0, 

zr  20 

a < r < » 

(33b) 

The  boundary  conditions 

in  polar  form  are 

u - u ■ A 
r 9 

0 < r < c 

(34a) 

a + u«  - 0 
r 9 

0 < r < d 

(34b) 

6 -a  -f(r) 
zr  z6 

c < r < 4 

(34c) 

CT  + O’  * 0 
zr  z6 

a < r < OB 

(34d) 

O +0  - 0 

zr  20 

6 < r < OB 

(34e) 

Next,  the  stresses  and  displacements  are  defined  In  the  form  of  the 
potentials  from  equations  (17),  (18),  and  (21). 


Equations  for  the  polar  displacements  and  stress  field  due  to  0^ 

and  0 are  given  in  [3].  Those  for  0 are  derived  by  coordinate  trans- 
2 3 

formation.  As  previously  noted,  n > 1 so  that 
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02  • „ ® Jl(5r)d5  cos  e 


T -S'/Aj 

« - j Aj(5)  t ■’  Jj(5r)d?  sin  e 


Making  uss  of  Chs  aquations  given  in  [3]  and  noting  that 


(35a) 


(35b) 


(35c) 
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•5*// ^ 


{^)  A^CS)  « 2^^(?r)j  d5  cos  e (36c) 


and,  with  some  regrouping  after  differentiation 


-52//V, 
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-5*//v, 
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(37f) 


The  above  stresses  and  dlsplacenents  are  now  used  with  the  boundary 
conditions  to  solve  for  the  and  subsequently,  the  stress  field. 

(b)  Boundary  Value  Problem  Solution 

From  the  form  of  the  boundary  conditions  (34)  it  becaoes  necessary 

to  define  (u  - iS  ),  (u  + fi),  (ct  - a )»  and  (a  + ) in  terms  of 

r 8 7 6 zr  z8  zr  z8 

the  A^(5). 

Conqtarlng  (26)  and  (36)  after  differentiation  it  is  seen  that 
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HAking  use  of  the  identity  [4]  that 
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it  follows  that 
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16 


For  (,$  ~ 9 ) and  (a  4-  ^ ) coimaring  (29)  and  (37)  after  differentiation 
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It  can  be  noted  at  this  point  that  from  syimietry  conaideratlona 

a - 0 on  the  entire  half 'space  surface,  and  a relationship  exists 
zz 

between  Aj^(§)  and  A2(^)>  sla^lifFlng  the  above  expressions.  From  (37f), 
on  z ■ 0,  equating  a to  zero  gives  rise  to 
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This  simplifies  the  expressions  at  hand  to 

\ 

1 

«• 

OB 

«r  + “q  “ J ^(5)  + A3(5)]  ^J^CIr)  d| 

0 

(44a) 

: 

: 

•% 

«B 

"r  ‘ “e  “ J 5J^(5r)  d5 

(44b) 

0 

i . 

I 

I 

Jk 

'rz  * ^ze  ■ '44  I [(^  * •=  ^ *3(1)] 

0*2  3 

■ i 

: 
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Define  a-l  + K,  b-c,,  f =-  + K ^ 1 , c ■ 
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Then 
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where  C(5)  - - | A^(|)  + 4^(5)  . D(§)  - - | A^(5)  - 43(5) 

and  d - ^ . e - %(d’4l),  f - %(d-^-l). 


The  boundary  conditions  (34)  can  now  be  written 


J CeC(?)  + fD(5)]  d|  - 0 


0 < r < c 


(49a) 
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(49b) 
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c J D(|)  5 J^(5r)  d5  - f(r) 


6 < r < a 


(49c) 
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C J D(5)  5 J^(?r)  d5  - 0 


S < r < 


(49d) 
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- c J C(§)  d5  - 0 


c<  r < 


(49e) 


Motivated  by  the  work  of  Westmann  [5]  and  Goodman  and  Keer  [6]  a 
solution  Is  assumed  In  the  form 


5 C(5)  - 0 


(50a) 


? D(5)  - J X^(t)  dt 


CL 

- 2f'  5^  J t\^(t)  J_3^(?t)  dt 
a 


(50b) 


where  f/  Is  the  coefficient  of  limiting  friction  to  be  subsequently 


Introduced. 


From  (50a)  It  Is  seen  that  (49e)  is  satisfied  Identically.  Sub- 
stituting (50b)  Into  (49d)  leads  to 


SjV')  ‘"=-“ 


which  Is  satisfied  fuitomatlcally. 

Substituting  (50b)  Into  (49c)  leads  to 


■fz  r ^ I dt. 


6 < r < fi 


At  this  point  the  nature  of  f(r)  is  examined.  As  noted  in  [2]  the 
pressure  distribution  under  the  area  of  contact,  0 < r < for  normal 


loading,  F , has  been  determined  as 
z 


p(r)  - — (a  - 


(53 


^ere  p , the  pressure  under  the  center  of  the  contact  region,  is 
o , . 


3P 


2ni 


(5^ 


Upon  application  of  a tangential  force,  P , sufficient  to  create 

X 


a slip  region,  c < r < S,  the  surface  traction  in  the  slip  region  may 
be  expressed  as 


<7_  - f ' p(r) 


£ < r < S 


(5f 


Then  by  (23)  and  the  fact  that  £(r)  * 2g(r), 


f(r)  - — ^ 


(5€ 


and  (51)  may  be  written 


B(t 


- i 4-  *r 


^ , r 0 (t)  Ct(t2-r^)’^]dt,  £ < r < fi 

dr  J o 


(57 
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PQ 

~ . KOotrlng  that  A (t)  must  satisfy  the  form  (57),  by 

Z fi  c o 

examining  p.  60  of  [4]  an  e3q>re8sion  for  0 (t)  is  found  to  be 

o 


0 (t)  - G(t^  - a^) 
o 


where,  defining  c by  (45), 


a. N ’‘j. 

^ 2 t \2c44>' 


Substituting  equations  (50)  into  (49b),  and  dividing  through  by  e 


gives  the  form 


-JT  4 . 1 dt 

V 2 e J (^2.^2)^ 


0 < r < c 


Examining  [4],  p.  67,  it  becomes  apparent  that  if  x^(t)  is  constant, 
(60)  will  be  satisfied.  With  Xj^(t)  ■ H,  integrating , 

■ VT  e “ »<^/2)  (61) 


and  H 


A ri 

e ''■/  TT 


. Therefore 


X,(.)  --tJJ 


satisfies  (49b).  It  is  easily  verified  that  (67a)  is  automatically 
satisfied  by  x^(^)  defined  above.  The  boundary  value 

problem  is  thereby  satisfied  by 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

1 

I 

I 

I 
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Xi(t)  - - S i/f 


(63) 


utilizing  the  identity  J , (5t)  - / -2-  co8(5t)  and  integrating  (50b), 
D(^)  is  found  to  be 


D(!)  - ^ I)  - it'eJJ  [m  . 2 iMISi 


. 2c  + 2 + (*2^2)  »ln^£S)j  (54) 


From  (47)  it  is  seen  that 


A^(5) 


^D(?) 


(65a) 


43(5)  - -%  D<5) 


(65b) 


and  by  (43) 


- i 


(65c) 


(c)  Auxiliary  Relationships 

Before  proceeding  to  the  evaluation  of  the  stresses  it  is  of  interest 
to  examine  two  conditions  key  to  the  determination  of  c,  the  radius  of 


the  no- slip  region,  and  the  relationship  between  P , the  tangentially 

X 


applied  force  and  a,  the  displacement  in  the  x-direction  on  the  surface 
under  the  area  of  contact. 


I.  - 


LI 
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In  the  work  of  Goodaen  and  Keer  [6],  upon  tangential  loading  and 

creation  of  the  slip  region,  £ < r < i,  a singularity  arises  In  a ^Ich 

zx 

can  be  eliminated  by  an  appropriate  choice  of  c.  To  examine  that 

phenomenon  In  this  analysis  attention  Is  turned  to  the  expresslcm  for 

a given  In  (32a).  Using  (46c)  and  (50a)  it  Is  seen  that  the  sum  of 
zx 

the  Fourier  coefficients  vanishes.  Therefore  any  singularity  In  a 

zx 

is  reflected  In  (a  ~ d ). 

zr  z0 

Performing  the  Indicated  Integration  and  differentiation  It  Is 
found  that 


zr  z6 


m ^ 2f 'G^C2(fi2.r2)%  - (fi^-r^)^] 


- Zf'G^^  (fi^-r2)(£2.r2)-% 


0 < t < £ 


As  the  radius  r goes  to  £ the  singularity  generated  In  (66)  can  be 


eliminated  If 


A - -2f'Gey^(fi^-£^) 


An  expression  for  the  radius  of  contact  under  normal  load,  p , 

z 

has  been  determined  (slmlllar  to  that  of  [2])  as 
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where  R is  the  radius  of  the  sphere,  and  6^  and  6^  are  quantities  to 


be  defined  later  \dilch  are  dependent  only  on  the  elastic  constants  of 
the  half-space.  Therefore  with  (69)  to  determine  R,,  (68)  can  be  solved 
for  c If Is  known. 

In  order  to  determine  the  relationship  between  the  displacement, 
ait  and  the  tangential  force,  P , (and  therefore  between  A and  F ) 


X “X 

equilibrium  must  be  considered.  On  the  surface  z « 0,  summing  forces 


In  the  x-dlrectlon  gives  that 


a 2rr 


-P 


J J 


<T  (r,0,o)  r d 0 d r 
zx 


(70) 


o o 


Again  using  (32a)  to  define  a on  the  surface,  with  the  sum  of 

ZX 


the  Fourier  stress  coefficients  equal  to  zero  as  noted  previously, 
performing  the  necessary  differentiation  and  Integration  on  the  dif- 
ference of  the  Fourier  stress  coefficients  it  is  found  that 


-Px  - 2c  j^-  A £ - f 


»%2  A 

a c - 


i-  ■ 1^)] 


(71) 


By  means  of  (68)  and  (71),  A (and  therefore  a)  and  c can  be 


calculated  given  i and  P^. 


1 


3.  Stress  Expressions 


(e)  Prellailnsry  Integration 


Before  proceeding  It  Is  necessary  to  connent  on  notation.  Since 
the  stresses  from  (37)  are  for  the  tangentially  loaded  spheres  and  will 
be  superposed  with  the  solution  of  Dahan  and  Zarka  [2]  for  the  normally 
loaded  case,  the  notation  to  follow  will,  as  much  as  possible,  be  con- 
sistent with  [2].  In  the  present  analysis  for  tangential  loading,  the 
elastic  constants,  c^^,  are  the  moduli  for  a given  transversely  Isotropic 
metal.  In  [2]  analysis  Is  carried  out  using  the  compliance,  a^^,  as  the 
elastic  constant  in  the  basic  equations.  The  relationship  between  the 
two  Is  herewith  stated  as: 


*11*33'*13 


*13  '*12*33 


*13*12"*11*13 


-2.2 

*11  *12 


’ "44 


— . (72) 


Further  relations  are  given  In  Appendix  A. 

It  Is  apparent  from  equations  (37)  that  several  key  Integrals  will 
be  Involved  In  the  stress  expressions.  From  equations  (65)  by  letting 


S • «3  - • 


It  can  b€  8««n  that  A^(;)  ■ 0(S)  vhare  1 ■ 1,2  or  3.  All  Beasal 

functions  In  (37)  are  In  sona  form  of  J^(Sr),  so  that  the  key  Integrals 


J D(5)e  V^(5r)  d5  , 


J 5 D(5)e  ^J^(5r)  d§  , 


^ 2 -5*//v. 

J 5 D(5)e  d§  , 


(75a) 


(75b) 


(75c) 


By  (64) 


J 5^D(§)e"^*^'^'^ij^<§r)  d| 


(•!9l 


-5*//v. 

8in(65)e  d§ 


-2f 'o/J  \-2i  J - co8ifi5i.J^-5*//vi 

+26  T raln.(i£)  . cosi65i1,-5*//Vi  J 
J L g,2  5 J 1 


2 2 - -?*//v,  . 

+ (6^-6^)  r sin(65)e  J^(5r)  d?} 
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j 5 D(5)  J^(5r)«  ‘d5 

o 

-2f 'oyl  [-2to^j(*>  + 2CD^^(C)  + (i2-«2)(.  ^ (o)] 

(80) 

B7  the  same  manner 

J 5 0(5)  J^(5r)e  ^d| 

o 

<» 

' (■  *4)  *1.1  [-“  “u<®  + **  ®u<*> 

+ (1^-6^)  J (£)]  (81) 

and 

r ”5*/yv£  / A 2N  / 1 \ 

J 0(5)  J^(5r)e  ^d5  . (-  41)  £ (o^^(£)  ^ x^^J(£)) 

-2f 'cyr  [-2fi  D^^(a)  + 26  D^^(£) 

+ (d2-fi2)  £ 1 )J 


29 


To  ovaluate  eh«  intagrsls  dofin^l  in  (77)  and  (78),  from  [2]  let 

2 2 

4 1 o *1  *»<»<»  *.* 


Y ^ (r^-  + a ^*^)^  + 4 


- %(r^  - A2  + + A^  y^) 


1 cr 


A-p^ 

S (A)  ■ axctaaf 


1 

\.i  <'•>  ■ ^ 


1 

I ‘rt) . 

1,1  r 


C^^(a)  - 1 - i (r  «)  + ^(«)) 


<=13<® 


- -L  (&^ 


Sa^zr 


^ (A  - ^ + s S‘)  S (A)  + s ^ (A) 

2A  2 1 1,1  4A  Ifl 


(83a) 


(83b) 


(83c) 


(84a) 


(84b) 


(84c) 


(85a) 


(85b) 
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1 / 1 1 \ 
i ‘ 1 - ® T . (fi))  (85c) 

12  2fi  ^ i l-.l  J-»l  J 


1/2  , «!  * \ 1 8l*  1 

“l3<*  ■ M ■ ' * — ) "l,l  * -5-  ^I.l 


Sl.l<«)+r/3 


D (£)  Is  not  given  In  [2]  and  Is  evaluated  by  definition  (78b)  and 
14 

integration  by  parts.  The  result  is 

Di4(^)  “ I J — Jl(|r)e  ^“^d? 

o ^ 


(85d) 


+ £ ] riln^  . cos^l  , 

^ O ^ -1  o 

- ¥ ] (86.) 


or,  by  the  previous  definitions 


D (ifi)  « ~ (^D  (A)  + T ^ (m)  ) + f C (m)  - --=•  D (m) 
14  4 \ 12  1,1  J 4 13  4 13 


(b)  Stresses  In  the  Half-Space 


(1)  Static  Case 

By  differentiation  and  substitution  of  the  appropriate  Integrals 
(80) -(82)  Into  (37)  the  stresses  In  the  half -space  become 


(86b) 


32 


T - {(-  H)(-  S [-^ 

+ 2«  Djj(«)  t (*V)  (-^ij.J(«))]} 


* (-^)  h {(-  ^ f)  ^ (\2‘«  + - “VI  • 

C-2S  D^^(a)  + 26  D^^(£)  + (S^-£^)£(D^^(fi)  + \^i(c))]} 


* '•i  {(■  H) « (“la'f^)  - “'=vT  • 

C-2S  (S)  +2£  dJ(£)  + (fi^-c^)£(p^^'(c)  + '(5))]} 


(87b) 


re 

sin  e 


c -c 

L 

' 2 / 3 


i("  ^ ' 4r  T ^(c))  -2f 'g  |’-2fi  D Ca) 

l\  e tt;\  d?  3,1  / ^ TT  L 32 


+ '■i  {(-  ^ ^ (“la*'^  * ^1 

+ 26  + (fi^-62)e 


+ (_12_U)  I,^  ir  e -a'cjji-u  D^^'(S) 


ZB  D^^(c)  + (a2.£2)£  ^ ^^^(fi))]} 


(87c) 


"^ze  C^(l+k4) 

iI5T  — ; — ^ 


‘j  ? {(-  f f>j.i(«  - “-'<=,/^ [-“  »J3(«  * 2e  »J3(«) 


- L if-  A 1'':  S.^/('£'>  - 2f'C  /T 


/V3  3 U e 


S3.1<^>  ’ 2f 'Gy  ^ C-2S  D33'(fi)  + 2fi  D '(S) 


+ (S^-c^)  s (c):^. 

i j 


(87d) 


^rz 

cos  9 /v. 


‘j{(-  *f)  “l3<«  '>33<‘) 


^ “33<*>  ♦ “ ”33<'> 


+ (a^-c^) 


(87e) 
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cos  6 


( Vj  ’ 


^ {(■  ?n)(-  s 


+ 2c-  D «)  + (4"-c-2)  (-  ^ } 


(87f) 


Equations  (87)  represent  e3q>resslons  for  the  stress  components  due 

to  a tangential  loading,  P^,  insufficient  to  cause  overall  sliding  of 

the  spheres.  The  results  of  Dahan  and  Zarka  [2]  for  the  stresses  due 

to  the  normal  load,  P , are 

z 


"rr  ■ '“o  .?Vi'  ■ “*2'*>) 


(88a) 


06  o L(s  -8^)(a'c  '-d')  V 1 2 12  2 1 22  J 


(=i'>2“i3<®>  - •2'‘l‘’23<®)’ 


(88b) 


^zz 


(88c) 


Q Mi  A 

rz  o (a^-a^yd 


(88d) 
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With  equations  (88)  superposed  with  equations  (87)  the  solution 
fox  the  stress  field  of  the  static  case  is  complete. 


(11)  Stresses  for  the  Sliding  Spheres 

From  the  static  case,  if  the  force  P is  gradually  increased,  the 

X 

size  of  the  no-sllp  region,  0 < r < c,  decreases  until  P reaches  a 

X 

value  of  f/p  , at  which  point  c becomes  zero  and  the  spheres  slide  over 
z 

the  entire  surface  of  the  contact.  The  effect  of  c-»0  is  analyzed 
and  generates  the  stress  expressions  in  the  half-space  due  to  slld- 


As  c-*0  it  can  be  seen  that 


SO  that  the  conqilete  stress  expressions  for  sliding,  including  con- 
tributions from  both  normal  and  tangential  loading,  can  be  written 


'11/-J-J2 


■) 


(89a) 
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L(—  * '12)  + {—;r-  ) h 


c,--c. 


+ (!il^)  (4f'G/T*)  CO.  e 


■ (“l-2‘>13<*>  - “2'>i”23'*>}| 


(89b) 


"re  * [-^—  '■3'>32<*>  ♦ 


(^)  \ 


”l4‘*> 


(!i2^)  L^Dj(g)]  (4f'G,/l  s)  .ic  e 


(89c) 


7>J3<*>  - ^ S^33'<*)]  («'<=./!  “)  e 


(89d) 


"r.  ■ [-7'*^’  "j“j3“>  (?)  V33<*>]  *)  S 

j ''  3 


Di2(a)-D22(a) 

’o  (8^-82)/d' 


(89e) 


t CO.  e 


+ p 


(•i-»j) 


(I 


where,  as  before,  J 

L3--M3. 


1,2  suomed,  i 


1,2,3  sunned, 


(ill)  Stresses  as  r-*0 

^ substituting  the  leading  terms  of  the  series  expansions  of  the 

Bessel  functions  Into  equations  (37)  the  stresses  can  be  examined  as 

the  radius  r approaches  zero  for  tangential  loading.  The  analysis  for 

normal  loading  has  been  performed  In  [2].  For  the  tangential  loading 

7 , a >a  , and  9 vanish  as  r-^0.  In  the  normal  case  9 vanishes 

rr  86  zz  rQ  rz 

and  9 and  9 do  not  arise  under  normal  loading. 
r8  ze 

For  the  case  where  tangential  loading  is  not  producing  sliding  of 
the  spheres,  the  superposition  of  the  stresses  for  r-»0  Is 


r. 


I 

T 
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®r8  " ° 


(90b) 


C44(I-rtc^)  ^ 2\/  S 


*0  Vv. 


{(-  i f)(^17)  -“'“71  [-“(^  - ‘-T 


(90c) 


-c^(l4k,) 


rz  Vv, 


“i  {(•  •^"'“71 


* ■ T ferv)]}  “•  ® 


(90d) 


■ "o  {-'  • n!  ■ t ■ “"“(A)] } 

12  1 2 


(90e) 


f^cre  1 * 1>2,3  suomed  and  k « 0 . 

3 

For  the  case  of  the  sliding  spheres  c-»0  and  equations  (90a,  b,e) 
remain  the  same.  Equations  (90c, d)  become 


a 

rz 


-C44(l-rts^) 


M. 


Si,i(fi))  COS  e 


(91a) 
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/v. 


(91b) 
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with  th«  same  aaaniatlon  convention  as  (90). 


(c)  Stresses  on  the  Surface 

On  the  surface  z « 0 the  exponential  term  In  (37)  goes  to  unity. 

The  Integral  expressions  (75)  are  likewise  simplified  and  evaluated. 

Performing  the  necessary  differentiation  and  substitution  In  (37), 
the  stress  egressions  on  the  surface  for  tangential  loading  are  obtained. 
The  stresses  on  the  surface  for  the  normal  loading  are  given  In  [2].  The 
superposed  results  give  the  stress  state  for  the  static  case  on  z - 0 as 
on 

0 < r < £ 


(92a) 


(92b) 


c ■ 0 (92c) 

to 

'r.  ‘ I {(■  • w) 

+ cm  e (924) 


^ . 


- I {-2f 'G JJ  C-2(fi^-r^)^]}  cos  Q 


- - ^(‘  - S)' 


and  on  d < r < 


.c,,k 


{(- - « V 


- t .ln'\D  - + r + (^ 


cr^-*V']}co,e.(21^X{(.Ai)[a 


c,,-c. 


.A 


it/ 


_2_^2,3/2  ,^2 


fi  ^ A-».  ^ O 


X.  2 .1 


44 


- i-^-h  {(-  ^ 

^ (f  ■ S ) 1 


<j  ■ a . - <j  « 0 
r*  z@  *2 


(94c) 


(94d) 


where  1 « 1,2,3  sumned,  j « 1,2  simnaed,  L » M , L ■ M . L •-M 

1 12  2*3*3 

before,  and  c Is  defined  by  (45). 

The  above  results  are  easily  speclaliaed  for  sliding 
as  follows: 


0 < r < d 

- [(^  - O - (^)  ^ cc,  e 


(95a) 


ee  • [(^  - =u)  (^'“71 " ')  - (^)  H s 


(95b) 
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4.  Numerical  Results 

In  order  to  show  the  effect  of  this  kind  of  anisotropy  on  the  stress 
state,  the  stresses  In  the  half-space  and  on  the  contact  were  examined 
numerically  for  magnesium  and  cadmium,  two  metals  studied  In  [2].  For 
these  calculations  the  coefficients  of  limiting  friction  f/.  >.  »25  and 
f/  ■ .50  were  used.  As  in  [1]  and  [2],  von  Mlses'  criteria  for  plastic 
yielding. 


{< 


*y  xz 


-O’  )^+(o  “O  )^+(a  “O  ) 

yy  Z2  ' zz  zx:  '•  xx  yy 


(97) 


was  calculated  and  lines  of  constant  J/p^  plotted  in  the  sphere,  on  the 
y a 0 plane,  and  cm  the  z » 0 surface  Inside  the  contact  zone.  Figures  4.1 
4,4  illustrate  the  results  in  the  half-space,  while  figures  4. 5-4. 8 illus- 
trate the  results  on  the  surface  for  r s;  S . Figure  4.9  represents  the 
static  case  on  the  surface  for  magnesium  with  f'  a .5,  and  £ > .5a,  and 
Is  Included  for  contrast  with  the  sliding  cases. 

It  Is  noted  In  [2]  that  It  Is  not  generally  possible  to  use  von 
Mlses'  criteria  to  predict  yielding  In  anisotropic  metals.  It  Is  plotted 
herein  as  a means  of  contrasting  the  stress  disturbance  In  the  trans- 
versely isotropic  metals  to  that  In  Isotropic  metals  as  displayed  In 
^1].  It  can  be  seen  that  magnesium  behaves  In  much  the  .same  manner  as 
the  Isotropic  case  In  [1],  while  cadmium's  elastic  properties  create  a 
more  unusual  situation.  Except  for  the  case  of  magnesium  with  f ' - .5 
(which  parallels  the  Isotropic  case  of  f ' • .5)  the  maximum  value  of 
J/p^  Is  found  on  z > 0 at  the  leading  edge  of  the  contact,  as  is  the  case 

In  [1]. 
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The  elastic  moduli  of  magnesium  and  cadmium  used  in  the  nusierical 
analysis  are  given  in  Table  4k  1.  Also  listed  are  the  moduli  of  steel 


for  V ■ .3 

Illustrating  values  for 

an  Isotropic  case. 

All  units  are 

N/m^. 

TABI£  4el 

ELASTIC  MODULI  FOR  MACIlESimf, 

CADMIUM, 

, AND  STEEL 

- 

Steel 

Magnesium 

Cadmium 

^®ir®33*‘^12"®13’ 

c ■%Cc  -c  ) 
44  ■ 11  12^^ 

c 

11 

5.857 

X 10 

1.092  X 10^^ 

2.691  X 10^^ 

c 

12 

2.501 

X 10 

3.976  X 10^° 

1.153  X 10^^ 

*^13 

2.079 

X 10 

10 

3.754  X 10 

11 

1.153  X 10 

®33 

6.110 

X 10 

4.602  X 10^*^ 

11 

2.691  X 10 

®44 

1.658 

10 

X 10 

10 

1.562  X 10 

10 

7.690  X 10 

1.678 

10 

X 10 

10 

3.472  X 10 

7.690  X 10^° 
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